Coulomb Effects in Spectral Density and Transverse Conductivity of Layered Metals 
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Renormalization of the Coulomb interaction in layered metals results in a strongly anisotropic 
plasma mode with low frequencies for small components of wave vector in the in-plane direction. 
Interaction of electrons with this mode was found to lead to an incoherent contribution to the electron 
spectral density which spreads up to large energies. In the superconducting state this reproduces the 
peak-dip-hump feature similar to that observed in layered high-T c superconductors. The incoherent 
part of the spectral density and plasmon-assisted electron transitions provide mechanisms for nearly 
linear conductivity in the stack direction at large voltages or frequencies. 

PACS numbers: 73.21.-b, 73.63.-b, 74.72.-h 
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Layered crystals can be considered as very anisotropic 
crystals with small interlayer transfer integral determin- 
ing the electronic bandwidth in the transverse direction. 
When such a crystal is pure enough so that its periodic- 
ity is not strongly disturbed the interlayer electron tran- 
sitions are expected to be coherent, i. e. the in-plane 
component of the momentum is conserved, in contrast to 
rough artificial tunnel junctions where it is not conserved. 
As it is well-known from the textbooks, electrons oscil- 
late in an electric field applied to an ideal crystal, and do 
not produce the dissipative current. A finite resistivity 
(and the real part of conductivity) in the stack direction 
appears due to electron scattering. At frequencies ex- 
ceeding the in-layer scattering rate, 1/r, the scattering 
becomes ineffective and conductivity is expected to de- 
crease as frequency increases. Similarly, the time needed 
for an electron to oscillate in the electric field becomes 
smaller than r if the voltage drop per a single layer be- 
comes larger than h/er, and the real part of the conduc- 
tivity must decrease with voltage increasing. However, 
one of the most studied layered metals, high-T c super- 
conductors, exhibit finite conductivity with nearly linear 
IV curves even at voltages much larger than the super- 
conducting gap d^ll, which in turn exceeds H/t. The 
real, dissipative part of the conductivity does not disap- 
pear at large frequencies as well Q . Similar behavior was 
observed in layered metal 2H-TaSe2 ||. 

Furthermore, the spectral density for electrons in high- 
Tc superconductors measured by means of angle-resolved 
photoemission spectroscopy contains the peak-dip-hump 
structure with an incoherent part spreading up to ener- 
gies far from the Fermi surface . The similar behaviour 
was observed also by means of tunneling spectroscopy 
f|. The presence of the incoherent part may account 
for the finite dissipative current at large frequencies 
As it was shown by Norman et al. the features in 
the spectral density may be modeled by coupling of elec- 
trons to some non-dispersive bosonic mode, presumably, 
of electronic origin. An alternative explanation was pre- 
sented by Chubukov and Morr who argued that the 



unusual superconducting properties of cuprates could be 
explained by a strong interaction between electrons and 
overdamped spin fluctuations peaked at some momen- 
tum. 

We suggest another mechanism for the discussed prop- 
erties that is related to general properties of layered met- 
als and is not confined to specific properties of high-T c 
superconductors. We calculate the effect of Coulomb in- 
teraction on the electron spectral density and conduc- 
tivity at large voltages and frequencies, and find that 
discussed features of layered metals can be induced by a 
strongly anisotropic plasma mode interacting with elec- 
trons. Such a mode must be inherent to layered materials 
because frequency of plasma oscillations with wave vec- 
tor perpendicular to conducting layers is proportional to 
the square of the interlayer transfer integral and, hence, 
to the transverse conductivity, the latter being small 
in highly anisotropic layered materials. The Josephson 
plasma mode in high-T c superconductors Jl^,[l3| can be 
considered as a manifestation of such a mode in the 
superconducting state. It is evident that similar mode 
must exist also in the normal state at frequencies larger 
than the scattering rate and at wave lengths shorter than 
mean-free-path, i. e. in the limit where a metal behaves 
as an ideal conductor, its response to electric field being 
in many respects similar to that of a superconductor. 

We consider a layered metal with the lattice period 
s in the direction perpendicular to the metallic layers 
coupled by small transfer integral t±. The Hamiltonian 
of the system is 
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7~tc + 7~Lbcs 



where n a is creation operator of electron with momen- 
tum p in the plane in the conducting layer n and with 
spin (j, m is the in-plane effective mass, Tic describes 
Coulomb interaction and Ti. bcs is the part of the Hamil- 
tonian of BCS type which leads to singlet s- or d-wave 
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pairing. We consider a discrete model in which the po- 
tential of Coulomb interaction of two electrons in layers 
n and n! separated by the in-plane distance ry also has 

a discrete form Vc = e 2 / ^Jr 2 + (n — n') 2 

transform reads 



V c =^e 2 /{q 2 + ql), 



(2) 



where q\\ is the in-plane wave vector, q^ = 
(2/s) sin (q±s/2) and \q±\ < tt / s is the wave vector ob- 
tained from the discrete Fourier transformation with re- 
spect to layer numbers. 

We neglect scattering of the electrons by impurities 
and phonons since typical energies involved in our study 
are large in comparison to Ti/t. We also ignore the reg- 
ular Coulomb effects like renormalization of the effec- 
tive mass and broadening of the quasiparticle peak in 
the electronic spectral function due to finite life-time of 
the quasiparticlcs. Instead we concentrate on the effects 
induced by plasma mode and its interaction with elec- 
trons. 

We calculate first the polarization operator within the 
RPA approximation adopted for the anisotropic case and 
find the renormalized Coulomb interaction, Vc at zero 
temperature. In the normal state and in the dynamic 
limit lu qnvp it reads 



VcO,q) = 
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where Q 2 



2e 2 VFPF I 's = Ane 2 n/m is the in-plane 
plasma frequency and u> p = (At _\_s / Tiv f)Q p <C Q p is the 
plasma frequency for wave vectors normal to the layers. 
Poles of the renormalized potential describe the spectrum 
of the plasma mode. The frequency of the mode starts 
from ujp and in a wide range of frequencies between oj p 
and £l p almost linearly increases with qu, o; q « q\\£l p /qi_. 
In the calculations below the leading contribution to in- 
tegrals comes from large values of q± ~ s . Thus 
the mode can be characterized by typical velocity values 
c ~ (ks)vf, where k 2 = fl 2 /2vp, 1/n is the Thomas- 
Fermi screening length. To simplify our calculations we 
assume that this length is smaller than the spacing s, 
ks > 1. This relation holds, e. g., for high-T c materi- 
als in which k ~ 10 mil -1 and s ~ 1.5 nm. Note that 
c» up, and, hence, we can use dynamic limit in (J3j) - 

There is some analogy between interaction of electrons 
with photons in quantum electrodynamics and interac- 
tion of electrons with plasmons in our problem. But while 
in the former case interaction was described by the small 
dimensionless parameter e 2 /(Tic) — 1/137 in the latter 
case a corresponding parameter e 2 /(Tic) ~ (e 2 / s) / (Tifl p ) 
can be either small or not depending on the plasma fre- 
quency and, hence, on the density of free charge carriers. 



We consider first the case of a relatively good metal with 
large carrier density when this parameter is small. Then 
we discuss what happens at smaller carrier concentration. 

The renormalized Coulomb potential for the super- 
conducting case in the quasiclassic limit, uj,qvp -C A, 
coincides with the potential in the normal state in the 
static limit, Vc = 47re 2 /(<7 2 + q]_ + n 2 ). Thus in the 
low-frequency limit in the superconducting state there is 
no singularities due to the plasma mode in the Coulomb 
potential. They recover only at w, qvF 3> A when the su- 
perconductivity is not important. However, the Joseph- 
son plasma mode at low frequencies with the spectrum 
similar to (S) does exist and is manifested by poles of 
the Green's functions describing fluctuations of the su- 
perconducting momentum (which is related to the vector 
potential). These fluctuations also contribute to the ef- 
fects we study here, but we do not study them here be- 
cause their contribution was found to be by (ks) 2 times 
smaller, than that from the Coulomb effects. 

Now we calculate the self-energy of electrons due to 
renormalized Coulomb interaction which acquires prop- 
erties of Green's function of bosons. The leading con- 
tribution to the mass operator is related to the poles of 
Vc describing the plasmons. This contribution comes 
from region of large energies. The resulting mass oper- 
ator does not decrease up to very large energies of the 
order of £l p . Using the calculated mass operator we find 
the electronic Green's function in the normal state at 
energies much smaller than fl p 
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where £ p = p 2 /2m — £f and F- is a slowly varying func- 
tion equal to 1 at £ < £ks and to 1/2 at e 3> £ks. We 
keep the second small term in the denominator because 
it determines the imaginary part of the Green's function 
at e ^ Formally the integral for go logarythmically 
diverges at large qii . But it is finite if we take into ac- 
count that integration over q\\ is limited by the Brillouin 
zone or consider the exact behavior of w q at q\\ ~ k. 
Then one gets go ~ g In sqF or go ~ glaKS, respectively, 
more detailed expression for go being dependent on the 
energy structure of the metal at energies far from the 
Fermi surface. 

Expression (Q) implies that the spectral function for 
electrons, A = ^sG/ir, consists of the suppressed quasi- 
particle peak decreasing with decrease of the carrier den- 
sity and of a broad feature extending up to the in-layer 
plasma frequency. For e <^ £l p it has a form 



A 



S(e-^p) , gF 



(5) 
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where we neglected broadening of the quasiparticle peak. 
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In the superconducting state the expression (@) is valid 
at e 3> A (res), while at e <C A the second term in the 
spectral function is small because it is determined by 
poles of the Coulomb potential, and the latter are not 
present at energies smaller than A. Thus the spectral 
density has a dip at energies e < A(res). Thus there is 
a peak-dip-hump structure similar to that observed in 
high-T c superconductors in (0, it) direction, which cor- 
responds to the maximum of the superconducting gap. 
However our analysis within the simple model with the 
isotropic electronic spectrum in the conducting plane 
cannot give a quantitative description of the spectral den- 
sity including the anisotropy related to the d-wave sym- 
metry of the superconducting order parameter. 

Note that though formally equations (H||) are valid at 
T = calculation at finite temperature shows that they 
can be used at temperatures smaller than plasmon fre- 
quencies which give leading contribution to the integral 
in (H|), i. e. at kT <C hvp/s, £p, f2 p . So, practically, the 
expression for Green's function is valid at any tempera- 
ture. 

Let us discuss now what happens at larger values 
of the parameter g describing the interaction strength. 
In this case we must solve the Dyson equations using 
renormalized Green's function. At this step one must 
take into account the vertex correction T. Vertex dia- 
grams similar to mass operator contain large contribu- 
tions from large energies due to plasma-mode poles in 
the renormalized Coulomb potential. However our anal- 
ysis shows that the vertex corrections do not lead to 
qualitative changes and with the accuracy of constant 
factors equation (Q) survives if expressions for g and 
<7o are assumed to contain the bare plasma frequency. 
But the actual plasma frequency is renormalized, and, 
hence, depends on the carrier density in different way, 
ftp = fi p o/(l + go), — ► w q o/(l + 9o), where subscript 
is related to bare quantities. 

To calculate current in the stack direction at large volt- 
ages we use non-equilibrium Green's functions in Keldysh 
diagram technique. Current density between the layers 
n and n + 1 is expressed via the Green's functions which 
are off-diagonal with respect to the layer indices: 



2et ± lr i\1 
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r 21 \ 
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2tt 3 



(6) 



where upper indices are related to the time contour, and 
lower indices describe the layer number. 

Diagrams describing two basic mechanisms of the ori- 
gin of linear conductivity at high voltages are presented 
in Fig.l. Solid and wavy lines denote renormalized 
Green's function and Coulomb potential, respectively, 
and cross denotes t±. These diagrams give strict results 
in the second order with respect to interaction and de- 
scribe current at large g qualitatively. Diagrams with 
larger number of Coulomb lines either give small correc- 
tions proportional to powers of V/£l p or correspond to a 



a) 



b) 



-X- 



n " n+1 n * n+1 

FIG. 1. Diagrams for non-diagonal with respect to layer 
number Green's functions, which give main contribution into 
a current in the stack direction. 



renormalization of the vertex which does not change the 
results qualitatively. 

For the coherent electron tunneling at large voltage 
drops per conducting layer (eV ^> %/t in the normal 
state and V ^> A in superconducting state) this process 
contributes to the current only due to the renormaliza- 
tion of the Green's function that results in the non-zero 
spectral density spreading up to large energy. In this 
limit, assuming eV 3> huj p (but eV 3> hCl p ) we get the 
linear I-V curve with conductivity 
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where S p is the area of the two-dimensional Brillouin zone 
of metallic layers. Comparing equation (^) with the stan- 
dard expression for the linear conductivity, a = lu 2 t/4it, 
one can see that the role of an effective scattering time 



in ai is played by 



1 + 



2tt P J 



In superconductors 



with d-wave pairing diagram la contribute to the con- 
ductivity at V < A due to transitions of quasiparticles 
via the superconducting gap Q| , the role of the scatter- 
ing time being played by where A is the maximum 
value of the gap. 

The processes given by the diagram lb describes con- 
ductivity due to plasmon assisted electronic transitions. 
Its contribution to the current has a self-explaining struc- 
ture 



j= I dkdpdp'e lks4e3t2±LJci 



TT 3 qlV 2 

[n p (l - n p /)(l + iV q ) - (1 - n p )n p 'iV q ] 

[S ' £p' + V + c^ q ) - S (£p -&-V + lo^)} 



(8) 



where n p and 7V q are Fermi and Planck distribution func- 
tions for electrons and plasmons. Calculating the inte- 
grals we find at eV 3> hu> p , A linear conductivity equal 
to 
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(9) 



Under our assumptions the first contribution is larger, 
than the second one, <J\ja% ~ Ks. 

Note that similar expressions (0-^) for the conductivity 
is valid in the linear regime at frequencies > %uj p , A, 
since the linear response at frequency lo is described by 
similar diagrams with voltage V replaced by frequency 
to. 
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Now we discuss a relation of our results to the ex- 
perimental data. The most published data are related 
to high-T c superconductors. The latter possess specific 
features which are not taken into account in our simple 
model. For example, we do not take into account pos- 
sibility of Van Hove singularities in the vicinity of the 
Fermi surface, and angle dependence of t±. Therefore, it 
is difficult to compare our results with the experimental 
data directly. However, qualitative consistency with the 
data is present. Our model gives qualitative explanation 
of the peak-dip-hump structure in the spectral density 
on the microscopic basis and provides the mechanism for 
the high voltage or frequency conductivity in the stack 
direction. 

The value of the conductivity of F^S^CaC^Oa; mea- 
sured Jl],f|,U at F > A is few times larger than that 
at small voltages, V < A. The calculated high voltage 
conductivity contains an effective scattering rate of the 
order of the Fermi energy, while conductivity at V < A 
the role of the scattering rate is played by the supercon- 
ducting gap A. For the isotropic t± this would imply that 
the conductivity at high voltages is smaller than at low 
voltages. But if the interlayer transfer integral is domi- 
nated by the transfer via apex oxygen, as was assumed 
in Ref. |15|, then the relation between these conductiv- 
ities may become the opposite, because conductivity at 
small voltages is dominated by the quasiparticles near 
the nodes of the superconducting gap where t± is mini- 
mal while to the conductivity at V > A contribute 
quasiparticles with any momentum directions, the con- 
ductivity being dominated by directions with large val- 
ues of the transfer integral. Thus though we cannot give 
quantitative description of the conductivity within our 
simple model, the results do not contradict to the exper- 
imental data. 
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